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Abstract 

A  non-destructive  inverse  method  is  applied  to  predict  the  global  temperature  distribution  at  the  interface  between  the  carbon  plate  and  the 
membrane  electrode  assembly  of  a  proton  exchange  membrane  fuel  cell  (PEMFC)  based  on  the  measured  temperature  data  on  the  outer  surface 
of  the  end  plate.  A  direct  problem  solver  is  built  to  provide  temperature  solutions  for  various  physical  conditions  and  is  then  incorporated 
with  the  conjugate  gradient  method  to  develop  an  inverse  method  for  temperature  prediction.  Three  specified  functions  are  taken  as  exact 
temperature  distributions  at  the  interfaces  and  are  predicted  by  using  the  inverse  method.  The  validity  of  the  inverse  approach  is  demonstrated 
by  investigating  the  effects  of  the  number  and  locations  of  the  measurement  points  and  the  experimental  uncertainty  on  the  accuracy  of  the 
predicted  solutions.  In  addition,  the  possibility  of  internal  temperature  prediction  based  on  the  temperature  data  on  the  edges  of  the  carbon 
plate  is  also  studied.  Results  show  that  the  inner  surface  temperatures  of  the  carbon  plates  in  contact  with  the  MEA  is  possible  to  be  determined 
based  on  the  measured  temperature  data  on  the  outer  surfaces. 

©  2004  Elsevier  B.Y.  All  rights  reserved. 
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1.  Introduction 

Proton  exchange  membrane  fuel  cells  (PEMFCs)  have 
been  developed  to  provide  electricity  for  a  variety  of  sta¬ 
tionary  and  transportation  applications.  In  order  to  produce 
commercially  viable  units,  it  would  be  necessary  to  pursue 
further  research  in  the  areas  of  water  and  thermal  manage¬ 
ment  [1-4],  cell  stacking  processes  [5-7],  designs  for  bipolar 
plates  and  flow-fields  [8,9],  and  fuel  processing  [10,11]. 

PEMFC  modeling  can  help  advance  this  technological  de¬ 
velopment  by  facilitating  the  understanding  of  the  physico¬ 
chemical  and  heat  transfer  phenomena  occurring  in  a  fuel 
cell.  Furthermore,  these  modeling  studies  provide  a  valuable 
tool  for  the  optimization  of  the  structure  as  well  as  the  ge¬ 
ometry  and  the  operation  conditions  of  the  fuel  cell  stacks. 
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For  example,  Yuan  et  al.  [12]  simulated  a  fully  developed 
laminar  flow  and  heat  transfer  process  in  the  carbon  plate 
channel  with  rectangular  or  trapezoidal  cross  section  and  a 
permeable  porous  wall.  Kumar  and  Reddy  [13],  on  the  other 
hand,  focused  on  the  improvement  of  the  performance  of  the 
PEMFC  by  optimizing  the  channel  dimensions  and  shape  in 
the  flow-field  of  bipolar/end  plates.  Fately,  Um  and  Wang 
[14]  used  a  computational  fuel  cell  dynamics  (CFCD)  model 
to  study  electrochemical  kinetics,  current  density  distribu¬ 
tion,  fuel  and  oxidant  flow,  and  multi-component  transport  in 
a  PEM  fuel  cell  with  cathode  carbon  plate  having  straight  or 
interdigitated  air  channels. 

Since  the  electrical  energy  produced  by  a  fuel  cell  is  ac¬ 
companied  by  a  nearly  equal  amount  of  thermal  energy,  the 
generated  heat  is  required  to  keep  the  working  temperatures 
of  the  fuel  cells  within  safe  levels.  Without  proper  thermal 
management,  the  performance  of  the  fuel  cells  may  not  be 
maintained.  A  too  high  cell  temperature,  even  in  confined  ar¬ 
eas,  may  lead  to  membrane  dehydration,  and  a  too  low  cell 


116 


C.-H.  Cheng,  M.-H.  Chang  /  Journal  of  Power  Sources  139  (2005)  115-125 


Nomenclature 

a,i  (i  =  0,  . . 5)  undetermined  coefficients  of  tempera¬ 
ture  function 

C  heat  capacity  (kJ(kg  °C)_  1 ) 

h  heat  transfer  coefficient  (W  (m2  °C)_1) 

H  height  (m) 

j  specific  objective  function 

/  objective  function 

k  thermal  conductivity  (W  (m  ° C)- 1 ) 

L  length  of  fuel  cell  (m) 

m ,  n  number  of  temperature  measurement  points  on 

outer  surface  in  x-  and  y-direction 
NX,  NY,  NZ  numbers  of  grid  points  in  x-,  y-,  and  z- 
direction 

g*  internal  heat  source  (W  m~3) 

r  random  number  varied  between  —  1  and  1 

t  time  (s) 

T  temperature  of  fuel  cell  (°C) 

T  simulated  experimental  outer  surface  tempera¬ 

ture  data  (°C) 

ra  ambient  temperature  (°C) 

W  width  of  fuel  cell  (m) 

x,  y,  z  Cartesian  coordinates  (m) 

Greek  symbols 
ft  step  size 

8  error  norm  (°C) 

y  conjugate  gradient  coefficient 

tc  search  direction 

6  exact  temperature  function  at  MEA/carbon 

plate  interface  (°C) 
p  density  (kg  m_ 3 ) 

<y  uncertainty  of  temperature  measurement 

Subscripts 

C  carbon  plate 

Cu  copper  plate 

e  end  plate 

ex  exact  solution 

g  gasket 

i,  j ,  k  grid  point  indices 

/,  /  temperature  measurement  point  indices 

Superscript 
n  iteration  step 


temperature  may  result  in  water  condensation  or  even  flood¬ 
ing.  In  order  to  ensure  proper  thermal  management  for  the 
fuel  cell,  it  is  essential  to  predict  the  temperature  distribution 
inside  the  fuel  cell.  However,  it  is  the  surface  temperature  of 
the  fuel  cell  that  can  be  readily  measured,  and  the  tempera¬ 
ture  inside  the  fuel  cell  is  usually  difficult  to  determine.  Or 


alternatively,  to  measure  the  internal  temperature  of  a  fuel 
cell,  one  may  use  destructive  methods  in  which  a  number  of 
temperature  sensors  are  inserted  into  the  fuel  cell  directly. 
Unfortunately,  the  destructive  methods  actually  damage  the 
structure  of  the  fuel  cell  to  a  certain  extent  and  disturb  the 
flow  and  current  fields  inside  the  fuel  cell.  In  these  circum¬ 
stances,  a  non-destructive  method  which  can  be  applied  in 
determination  of  the  internal  temperature  of  a  fuel  cell  is  re¬ 
quired. 

In  this  study,  an  inverse  concept  is  proposed  for  obtaining 
the  global  temperature  distribution  inside  the  fuel  cell  at  the 
interface  between  the  ME  A  and  the  carbon  plate.  This  method 
features  the  major  advantage  of  producing  no  destructive  ef¬ 
fects  on  the  fuel  cell  structure.  With  the  present  method,  the 
interface  temperature  distribution  between  the  carbon  plate 
and  the  MEA  is  determined  based  on  the  temperature  mea¬ 
surement  of  the  outer  surface  of  the  end  plate.  The  validity  of 
the  present  inverse  method  is  tested  by  a  number  of  simulated 
cases. 

Fig.  1  shows  the  schematic  of  a  single  cell  of  the  PEMFC 
and  its  physical  model  as  well.  The  PEMFC  for  a  single  cell 
shown  in  this  figure  consists  of  a  polymer  electrolyte,  which 


Simulated  zone 


©  Membrane  electrode  assembly 
®  Flow  channels 
©  Carbon  plate 


©  Copper  current  collector 
©  Gasket 
®  End  plate 


.V 


Temperature  measurement  point  on  Grid  points  (NXxNY)  on  outer 
outer  surface  of  end  plate  ( w  xh)  surface  of  end  plate 


Fig.  1.  Physical  model  for  a  single-cell  PEMFC. 
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Table  1 


Fixed  dimensions,  material  properties,  and  surface  conditions  of  the  test 
cases 


Carbon  plate 

Copper  plate 

Gasket 

End  plate 

L 

0.21 

0.21 

0.21 

0.21 

W 

0.21 

0.21 

0.21 

0.21 

H 

0.003 

0.003 

0.0003 

0.004 

k 

95 

380 

0.17 

200 

h 

10a 

10a 

10a 

10b 

Ta 

25a 

25a 

25a 

25b 

a  Edge  convection. 
b  Outer  surface  and  edge  convection. 


behavior;  therefore,  the  three-dimensional  steady-state  heat 
conduction  equation  for  all  the  solid  materials  is  expressed 
as  follows: 


d1 2 3T  d2T  d2T 
dx2  3 y2  3 z2 


Note  that  the  MEA  is  not  necessarily  to  be  included  in  the 
present  modeling  since  the  solution  domain  is  confined  be¬ 
tween  the  outer  surface  and  the  ME  A/carbon  plate  interface. 
The  boundary  conditions  associated  with  Eq.  (2)  are: 

(a)  On  outer  surfaces: 


is  sandwiched  between  two  electrodes  and  two  gas  diffu¬ 
sion  layers  to  form  a  five-layer  MEA.  The  five-layer  MEA 
is  placed  between  two  carbon  plates  with  machined  grooves 
that  provide  flow  channels  for  distributing  fuel  and  oxidant 
gases.  The  inner  face  of  each  of  the  carbon  plates  is  in  contact 
with  the  MEA,  and  the  temperature  distribution  on  the  inner 
surface  of  the  carbon  plate  is  to  be  determined  in  the  study. 
The  outer  faces  of  the  carbon  plates  are  in  contact  with  two 
respective  copper  current  collectors,  which  are  insulated  by 
two  gaskets  from  the  outside.  Finally,  the  structure  is  clamped 
tightly  between  two  end  plates.  The  geometric  conditions  of 
the  materials  of  the  components  considered  in  the  study  are 
listed  in  Table  1 .  As  the  fuel  cell  is  in  operation,  heat  is  gener¬ 
ated  from  the  electrochemical  reaction  and  conducted  toward 
the  outer  surfaces,  and  then  dissipated  to  the  ambient  air  at 
temperature  Ta  by  convection.  The  heat  transfer  coefficient 
on  all  the  outer  surfaces  of  the  fuel  cell  is  denoted  by  h.  An 
array  of  temperature  sensors  is  installed  on  the  outer  surface 
of  the  end  plate  or  on  the  edge  surfaces  of  the  carbon  plate 
to  measure  the  surface  temperatures.  However,  note  that  the 
temperature  distribution  on  the  outer  surface  of  the  end  plate 
can  also  be  detected  by  using  an  infrared  radiation  thermal 
image  system.  The  effects  of  the  number  and  locations  of 
the  temperature  measurement  points  and  the  experimental 
uncertainty  are  among  the  major  issues  investigated  in  this 
study. 


±k^=h[T-Ta]  (3a) 

on 

where  n  is  the  coordinate  normal  to  the  individual  outer  sur¬ 
faces. 

(b)  At  interface  between  any  two  solid  layers: 


T\  =  T2 


(3b) 

(3c) 


where  the  subscripts  1  and  2  denote  any  two  solid  layers  in 
contact. 

Eqs.  (2)  and  (3)  can  be  transformed  into  a  set  of  si¬ 
multaneous  algebraic  equations  for  the  temperature  solu¬ 
tions  at  all  the  grid  points  by  using  the  finite-difference 
method.  With  the  help  of  the  successive-over-relaxation 
method  (SOR)  [15],  the  convergent  numerical  solution  for 
the  three-dimensional  temperature  distribution  within  the  en¬ 
tire  fuel  cell  (with  NX  x  NY  x  NZ  grid  points)  can  be  ob¬ 
tained.  Note  that  NZ  is  the  sum  of  NZc,  NZcu,  NZg,  and 
NZe,  the  numbers  of  grid  points  in  Z-direction  within  car¬ 
bon  plate,  copper  current  collector,  gasket,  and  end  plate. 
Fixed  physical  variables  used  in  computation  are  listed  in 
Table  E 


2.  Optimization  method 

2.7.  Heat  conduction  problem  solver 

The  materials  of  the  components,  including  the  carbon 
plate,  the  copper  plate,  the  gasket,  and  the  end  plate,  used  in 
the  fuel  cell  are  assumed  to  be  homogeneous  and  isotropic. 
Therefore,  heat  conduction  within  these  solid  layers  is  gov¬ 
erned  by  the  following  partial  differential  equation: 

dT  9 

pC —  =  kV2T  +  q  (1) 

where  cf  denotes  the  internal  heat  source;  p,  C,  and  k  are 
the  density,  heat  capacity,  and  thermal  conductivity  of  the 
solid  materials,  respectively;  and  T  is  the  temperature.  The 
present  approach  is  mainly  concerned  with  the  steady  thermal 


2.2.  Conjugate  gradient  method 

The  temperature  measurement  data  on  the  outer  surface  of 
the  end  plate  are  denoted  by  T j  j  (/  =  1  ~  m,  and  /  =  1  ~  n). 
Let  Ti ,;,nz  denote  the  numerical  temperature  solution  at  the 
grid  point  (i,j,  NZ)  with  i =  1  ~  NXJ  = 1 ~  NY,  and  k  =  NZ, 
which  is  on  the  outer  surface  of  the  end  plate.  The  number  of 
the  measured  temperature  points  ( m  x  n)  can  be  less  than  or 
equal  to  the  number  of  grid  points  (NX  x  NY)  on  the  outer 

Table  2 

Exact  solutions  at  MEA/carbon  plate  interface  for  the  test  cases 
CaseExact  solution 

1  6(x,  y)  =  80  +  lOx  —  80v2  +  20_y  —  80y2  —  5xy 

2  0(x,  y)  =  [sin(xy7r/L  W)]2  +  2  sin(x7i/L)  +1.5  sin(_y7r/  W)  +  81 

3  9{x,  y)  =  e-(*-0-12)2e-<?-°-1)2/{(0.012  +  O.l/)[cosh(xy;r/0.441)]0'4) 
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surface.  An  objective  function  of  the  present  inverse  problem 
is  thus  defined  as 

m  h 

7  =  EE  (T>-wz  -  TuY  (4) 

i(=I)  j(=J ) 


Note  that  the  temperatures  at  the  grid  points  on  the  outer 
surface  of  end  plate  (7}  7  nz)  are  numerically  provided  by 
solving  the  heat  conduction  equation,  and  the  temperatures 
at  the  measurement  points  (77?  /)  are  obtained  by  experiment. 
During  the  iterative  process,  the  temperature  distribution  at 


*  [m] 


x[m] 


x[m] 


(a)  Predicted  temperature 
distribution 


x[m] 

(b)  Measurement  point  arrangement 


Fig.  2.  Temperature  predictions  with  different  measurement  point  arrangements  on  outer  surface  of  the  end  plate,  for  case  3  at  o  =  0.1. 
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the  MEA/carbon  plate  interface  is  updated  iteratively  toward 
minimization  of  the  objective  function.  Minimization  of  the 
objective  function  means  minimization  of  the  difference  be¬ 
tween  Tt  j  nz  and  T j  j.  Therefore,  the  temperature  distribu¬ 
tion  at  the  MEA/carbon  plate  interface  is  obtained,  leading  to 
a  minimal  difference  between  the  numerical  and  experimen¬ 
tal  outer  surface  temperature  data. 

The  temperature  distribution  T(x,  y)  at  the  MEA/carbon 
plate  interface  is  approximated  by  a  polynomial  expression 
as 

9  9 

T(x,  y)  =  ao  +  a\x  +  a2X  +  a^y  +  a\y  +  a$xy  (5) 

where  ao,  a\,  a2,  <23,  <24,  and  a$  are  the  undetermined  coeffi¬ 
cients  to  be  optimized  in  the  iterative  process.  Different  com¬ 
binations  of  these  coefficients  represent  different  temperature 
distributions  among  which  the  optimal  interface  temperature 
distribution  may  be  found.  In  other  words,  in  the  optimization 
process,  the  coefficients  ao,  a\,  a2,  a3, 04,  and  a$  are  updated 
iteratively  toward  the  minimization  of  the  object  function. 

The  minimization  of  the  objective  function  J  is  accom¬ 
plished  by  using  the  conjugate  gradient  method.  The  con¬ 
jugate  gradient  method  evaluates  the  gradient  functions  of 
the  objective  function  and  sets  up  a  new  conjugate  direction 
for  the  updated  undetermined  coefficients  with  the  help  of 
a  direct  numerical  sensitivity  analysis.  The  method  of  sen¬ 
sitivity  analysis  by  direct  numerical  differentiation  was  first 
presented  by  Cheng  and  Wu  [16].  The  inverse  approach  has 
been  proved  efficient  in  the  applications  for  problems  with 
conduction  [16]  and  convection  [17].  More  recently,  it  has 
been  used  to  identify  the  profile  of  an  unknown  interior  ob¬ 
ject  in  a  solid  body  [18].  In  general,  the  convergence  can 
be  attained  in  a  finite  number  of  iterations.  The  approach  is 
described  briefly  in  the  following. 

The  following  is  a  summary  of  the  sequence  of  the  iterative 
optimization  process  for  searching  the  optimal  temperature 
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0 


'  I  ' 

0.05 


j  L 

0.1 


J  I  I  I  I  I  I  L 

0.15  0.2 


x  [m] 


of  the  MEA/carbon  plate  interface: 

(1)  Collect  the  experimental  temperature  data  on  the  outer 
surface  of  the  end  plate,  T j  j. 

(2)  Make  an  initial  guess  for  the  temperature  distribution 
at  MEA/carbon  plate  interface  by  giving  a  set  of  initial 
values  for  the  undetermined  coefficients  ao,  a\,  a2,  a3, 
a4,  and  a 5. 


/. 


Fig.  3.  Predicted  temperature  profiles  along  y  =  0.105  m  based  on  the  third 
and  the  fourth  arrangements  considered  in  Fig.  2. 


Fig.  4.  Effects  of  experimental  uncertainty  on  accuracy  of  temperature  pre 
diction  for  case  1 . 
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Table  3 


Obtained  coefficients  at  different  experimental  uncertainties  based  on  41  x  41  measurement  points  for  case  1 


Coefficient 

Experimental  uncertainty  (a) 

0 

0.1 

0.5 

1.0 

a0 

80.000 

80.021 

80.083 

79.892 

a\ 

10.000 

9.253 

9.219 

9.976 

02 

-80.000 

-76.850 

-76.172 

-73.074 

20.000 

20.170 

18.237 

23.491 

<24 

-80.000 

-81.406 

-72.575 

-91.338 

<25 

-5.000 

-3.760 

-3.578 

-22.687 

e[°C] 

4.028E— 6 

1.074E— 2 

2.624E— 2 

7.317E— 2 

(3)  Use  the  heat  conduction  problem  solver  to  obtain  the 
three-dimensional  temperature  distribution  within  the 
entire  fuel  cell,  and  read  7)  j  nz- 

(4)  Calculate  the  objective  function  /  by  Eq.  (4).  When  the 
objective  function  reaches  a  minimum,  the  solution  pro¬ 
cess  is  terminated.  Otherwise,  proceed  to  step  (5). 

(5)  Perform  the  direct  numerical  sensitivity  analysis  [16]  to 
determine  the  gradient  functions  dJn  /dai  (i  =  0,  1, 

5). 

(6)  Calculate  the  conjugate  gradient  coefficients  y1/  and  the 
search  directions  7r"+1  with 


(dJ/dcuT 


2 


(dj/dciif-1 


71  +  1 


dJ 


n 


da i 


+  yf+rf ,  i  =  0,  1, . . . ,  5 


(7) 


(7)  Calculate  the  step  sizes  /3;  (/  =  0,  1,  . . .,  5)  that  lead  to 
dJn+l /dpi  =  0,  by  solving  the  following  simultaneous 
equations: 


dTn  dTn 
01  Uh NZ  NZ 


da1/ 


dao 


n 


+ 


dTn  dTn 
01  Uh NZ  NZ 


da1/ 


da  i 


+ 


n 


+  @5  ^5 


X 


m  n 


dTn  dTn 
01  i,j, NZ  NZ 


7  =  1  7=1 


da1/ 


das 


m  n 

=  EE 

7  =  1  7=1  L 


7,7,NZ 


-  TU  ) 


dTn 

01  Uh  NZ 
da1/ 


i  =  0,  1,...,5 


(8) 


(8)  Update  the  undetermined  coefficients  with 


ai 


71+1 


=  a 


n 


a  perturbation  to  the  exact  temperature  distribution  on  the 
outer  surface.  That  is, 

T  I,J  =  (7+Nz)ex  +  <WjJ  (10) 

where  (7y;  Nz)ex  is  the  exact  temperature  distribution  at  the 
grid  points  located  at  the  same  positions  as  the  same  mea¬ 
surement  points  on  the  outer  surface  of  the  end  plate  with 
a  specified  interface  temperature  distribution  of  interest;  cj 
is  a  random  number  evenly  distributed  between  —1  and  1 
and  is  provided  by  a  random  number  generator;  and  a  is  a 
value  given  to  simulate  the  uncertainty  in  the  measurement. 
Note  that  at  a  =  0,  the  exact  outer  surface  temperature  solu¬ 
tion  (Tyy  Nz)ex>  is  used  directly  for  the  interface  temperature 
prediction. 

3.  Results  and  discussion 

In  order  to  investigate  the  validity  of  the  present  approach, 
three  different  cases  with  different  forms  of  the  exact  solution 
are  specified  to  be  predicted,  as  can  be  seen  in  Table  2.  There 
are  three  major  parameters  affecting  the  accuracy  of  temper¬ 
ature  prediction  considered  in  this  report.  They  are:  (1)  the 
arrangement  and  positions  of  the  measurement  points;  (2) 
the  experimental  uncertainty;  and  (3)  the  number  of  mea¬ 
surement  points.  Meanwhile,  the  possibility  of  temperature 
prediction  based  on  the  temperature  data  on  the  edge  surfaces 
of  the  carbon  plate  is  also  tested.  Results  obtained  by  using 
the  temperature  data  on  the  carbon  plate  edges  are  provided 
in  order  to  compare  with  the  data  obtained  by  using  the  tem¬ 
perature  data  on  the  outer  surface  of  the  end  plate.  Typically, 
41  x  41  x  35  grid  points  in  the  solution  domain  are  adopted 
after  a  grid-independence  test  has  been  performed.  That  is, 
NX  =  41  and  NY  =  41.  On  the  other  hand,  the  number  of 
measurement  points  on  the  outer  surfaces  of  the  end  plate  or 
on  the  carbon  plate  edges  is  varied  as  a  parameter. 

3.1.  Using  temperature  data  on  outer  surface  of  the  end 
plate 


and  go  to  step  (3). 

2.3.  Simulated  experimental  data 

In  this  study,  the  temperature  measurement  data  on  the 
outer  surface  of  the  end  plate  (T u)  are  simulated  by  adding 


There  are  four  possible  arrangements  for  the  installation 
of  the  measurement  points  on  the  outer  surface  of  the  end 
plate,  which  are  shown  in  Fig.  2.  The  arrangements  of  mea¬ 
surement  points  are  plotted  in  the  right  portion  of  this  fig¬ 
ure  and  the  corresponding  predicted  temperature  contours 
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are  shown  in  the  left  portion,  for  case  3  given  in  Table  2  at 
a  =  0.1.  In  case  3,  the  exact  temperature  function  is  spec¬ 
ified  by  (*—0-12)2£—(y— 0.1)2/{(0. 012+0.  Iy5)[cosh(xy7r/O.441)]04} . 

As  clearly  shown  in  this  figure,  with  the  first  and  the  sec¬ 
ond  measurement  point  arrangements,  the  approach  leads  to 
unrealistic  temperature  solutions  for  the  ME  A/carbon  plate 
interface.  Nevertheless,  the  predicted  solutions  of  the  third 


and  the  fourth  measurement  point  arrangements  appear  to  be 
identical  to  the  exact  solution.  In  order  to  investigate  the  dis¬ 
crepancy  between  the  solutions  obtained  from  the  third  and 
the  fourth  measurement  point  arrangements,  the  temperature 
profiles  along  y  =  0. 105  m  are  plotted  and  are  compared  with 
the  specified  exact  temperature  function  in  Fig.  3.  The  solid 
curve  in  this  figure  indicates  the  exact  temperature  solution.  It 
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Fig.  5.  Effects  of  experimental  uncertainty  on  the  accuracy  of  temperature  prediction  for  cases  2  and  3. 
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Fig.  6.  Numerical  data  of  best  fit  to  the  simulated  experimental  data  for 
temperature  distribution  on  outer  surface  of  end  plate  along  y  =  0.105  m  for 
cases  1,  2,  and  3  at  a  =  0.1. 

is  found  that  the  third  measurement  point  arrangement  results 
in  a  greater  error  in  prediction  than  the  fourth  arrangement 
since  the  former  has  much  less  temperature  measurement 
points  than  the  latter.  As  the  accuracy  of  the  temperature  pre¬ 
diction  is  desired,  a  uniform  arrangement  of  the  measurement 
points  like  the  fourth  arrangement  is  essential.  However,  if  the 
experiment  cost  is  of  major  concerns,  the  third  arrangement  is 
still  acceptable.  In  the  present  study,  the  fourth  arrangement 
is  regarded  as  standard  arrangement  for  measurement  points. 

Fig.  4  displays  the  effects  of  the  uncertainty  of  the  exper¬ 
imental  temperature  data  (a)  on  the  temperature  prediction. 
In  this  figure,  the  exact  temperature  solution  is  specified  with 
case  1  of  Table  2.  The  temperature  prediction  is  based  on 
a  uniform  distributed  41  x  41  measurement  points  on  the 
outer  surface  of  the  end  plate.  In  this  case,  the  exact  solution, 
0(x,y )  =  80  +  l(k  —  80x2  +  20y  —  80y2  —  5xy,  is  a  polynomial 
function.  In  the  upper  plot  of  this  figure,  the  exact  tempera¬ 
ture  function  of  case  1  is  shown.  The  predicted  temperature 
profiles  along  x  =  0.105  m  and  y  =  0.105  m  are  displayed  at 


11x11  6x6  3x3 


y  [m] 


Fig.  7.  Effects  of  number  of  temperature  measurement  points  on  outer  sur¬ 
face  of  end  plate  on  predicted  temperature  profile  along  x  =  0. 105  m  for  case 
1  at  o  =  0.1. 


a  =  0  and  1.0.  In  the  lower  two  plots  of  this  figure,  the  solid 
curves  indicate  the  exact  solutions  and  the  dashed  and  dash- 
dot  curves  represent  the  predicted  temperature  profiles.  It  is 
found  that  at  a  -  0,  an  accurate  temperature  prediction  can 
be  obtained.  However,  when  a  higher  uncertainty  of  o  -  1.0 
is  introduced  into  the  experimental  temperature  data,  the  de¬ 
viation  between  the  predicted  and  the  exact  temperature  data 
is  greatly  increased. 

To  quantify  the  accuracy  of  the  temperature  predictions 
for  the  ME  A/carbon  plate  interface,  an  error  norm  is  defined 
based  on  the  discrepancy  between  the  predicted  and  the  exact 
temperature  distributions  at  the  ME  A/carbon  plate  interface. 
That  is, 


s  = 


ij  \(Tij,\)  Oij | 

(NX  x  NY) 


(11) 


where  T(jt  i  and  Oij  denote  the  predicted  and  the  exact  tem¬ 
perature  functions  at  grid  (/,  j)  at  the  ME  A/carbon  plate  in¬ 
terface  and  (NX  x  NY)  are  the  number  of  grid  points  on  the 
interface  plane.  Table  3  shows  the  obtained  undetermined  co¬ 
efficients  and  the  error  norms  at  various  uncertainties.  In  this 
table,  the  uncertainty  is  varied  from  0  to  1.0.  It  is  observed 
that  the  error  norm  (s)  of  temperature  prediction  increases 
with  the  uncertainty,  and  the  predicted  coefficients  exhibit 
greater  error  as  o  is  elevated.  The  error  norm  s  reaches  ap¬ 
proximately  0.073  °C  at  a  -  1.0.  This  reveals  the  accuracy  of 
the  present  approach. 


Table  4 


Specific  objective  functions  at  different  experimental  uncertainties  for  cases  1,  2,  and  3 


Specific  objective  function 

Case 

Experimental  uncertainty  (a) 

0 

0.1 

0.5 

1.0 

J 

1 

1.535E— 11 

3.277E— 3 

8.065E— 2 

3.251E— 1 

2 

2.093E— 3 

5.357E— 3 

8.160E— 2 

3.360E— 1 

3 

1.466E— 4 

3.438E— 3 

8.356E— 2 

3.291E— 1 
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Table  5 


Obtained  coefficients  with  different  numbers  of  measurement  points  on  outer  surface  of  end  plate  for  case  1  at  a  =  0.1 


Coefficient 

Measurement  points  (m  x  n) 

21  x  21 

11  x  11 

6x6 

3x3 

a0 

79.981 

80.010 

79.977 

79.864 

a\ 

9.898 

9.601 

12.403 

9.609 

02 

-80.111 

-77.594 

-90.349 

-72.643 

a3 

20.531 

19.720 

18.137 

20.965 

04 

-83.164 

-78.053 

-69.510 

-73.430 

a5 

-3.502 

-5.521 

-8.276 

-19.134 

s[°C] 

1.023E— 2 

1.045E— 2 

4.181E— 2 

7.275E— 2 

Fig.  5  provides  the  results  of  predictions  for  the  non¬ 
polynomial  exact  temperature  funcions  with  cases  2  and 
3  based  on  41  x  41  measurement  points  at  different  un¬ 
certainty.  The  data  shown  are  similar  to  those  given  in 
Fig.  4.  However,  in  Fig.  5,  it  is  seen  that  even  at  a  -  0, 


for  case  2  and  3,  the  exact  temperature  curves  cannot  be 
fitted  exactly  by  the  polynomical  temperature  distribution 
like  case  1.  The  error  of  the  temperature  prediction  with 
case  2  is  greater  than  that  with  case  3.  The  results  imply 
that  the  magnitude  of  error  is  also  dependent  on  the  features 


T  [°C]  T  [°C] 


z 


B 

T  [°C]  T  [°C] 


Fig.  8.  Numerical  data  of  best  fit  to  the  simulated  experimental  data  for  temperature  distributions  on  edge  surfaces  of  carbon  plate  along  z  =  0.0015  m  for  case 
2  at  cr  =  0.1. 
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of  temperature  distribution  at  the  MEA/carbon  plate  inter¬ 
face. 

On  the  other  hand,  a  specific  objective  function  is  de¬ 
fined  in  order  to  evaluate  the  convergence  performance  of 
the  method.  The  specific  objective  function  is  defined  as 


J 

(m  x  n) 


(12) 


where  m  x  n  is  the  total  number  of  measurement  points  on 
the  measured  surface. 

Fig.  6  shows  the  comparison  between  the  optimal  numer¬ 
ical  data  and  the  simulated  experimental  data  for  the  temper¬ 
ature  distribution  on  the  outer  surface  of  the  end  plate  along 
the  y  =  0.105  m  for  cases  1,  2,  and  3  at  o  =  0.1.  From  the 
temperature  profiles  along  y  =  0.105  m  shown  in  this  figure, 
it  is  seen  that  a  best  fit  can  be  accurately  yielded  for  each  set 
of  the  simulated  experimental  data.  Table  4  gives  the  specific 
objective  function  at  various  experimental  uncertainties  for 
cases  1,  2,  and  3.  Based  on  the  data  listed  in  this  table,  it  can 
be  observed  that  an  increase  in  uncertainty  leads  to  a  higher 
j,  which  represents  a  worse  convergence. 

The  cost  of  the  experiment  is  reduced  when  the  number 
of  temperature  measurement  points  is  decreased.  Therefore, 
the  effects  of  the  number  of  the  measurement  points  on  the 
outer  surface  on  the  accuracy  of  the  temperature  prediction 
for  the  MEA/carbon  plate  interface  are  examined,  and  part 
of  the  results  are  presented  in  Fig.  7.  In  this  case,  the  consid¬ 
ered  number  of  measurement  points  in  the  x-  and  y-directions 
( m  x  n)  on  the  outer  surface  of  the  end  plate  is  assigned  to  be 
11  x  11  ,6  x  6,  and  3  x  3.  The  corresponding  positions  of  the 
measurement  points  are  plotted  in  the  figure.  Meanwhile,  the 
solid  curve  indicates  the  exact  temperature  profile  to  predict 
and  the  dashed,  dashdot,  and  dash-double-dot  curves  indi¬ 
cate  the  results  obtained  from  the  11  x  11, 6x6,  and  3  x 
3  measurement  points,  respectively.  It  is  observed  that  1 1  x 
1 1  measurement  points  lead  to  an  accurate  prediction  of  the 
MEA/carbon  plate  interface  temperature.  The  discrepancy  in 
the  predicted  temperature  distributions  between  the  exact  and 
the  predicted  curves  is  increased  as  the  number  of  measure¬ 
ment  points  is  reduced. 

The  dependence  of  the  temperature  predictions  on  the 
number  of  measurement  points  is  also  shown  numerically 
in  Table  5.  The  obtained  coefficients  of  the  interface  temper¬ 
ature  function  are  compared  with  the  exact  ones.  According 
to  the  results  provided  in  Table  5,  it  is  found  that  as  the  num¬ 
ber  of  measurement  points  is  decreased  from  21x21to3x 
3,  the  error  norm  s  is  elevated  from  0.01  to0.073°C. 


3.2.  Using  temperature  data  on  edge  surfaces  of  the 
carbon  plate 

The  possibility  of  internal  temperature  prediction  based 
on  the  temperature  data  on  the  edge  surfaces  of  the  carbon 
plate  is  also  studied.  Fig.  8  shows  comparison  between  the 
numerical  data  of  best  fit  and  the  simulated  experimental  data 
for  temperature  distribution  on  the  edge  surfaces  of  the  car- 
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Fig.  9.  Effects  of  number  of  temperature  measurement  points  on  edge  sur¬ 
faces  of  carbon  plate  on  predicted  temperature  profile  along  _y  =  0. 105  m  for 
case  3  at  a  =  0.1. 

bon  plate  along  z  =  0.0015  m  for  case  2  at  o  -  0.1.  In  this 
figure,  data  on  four  edges  (A-B,  B-C,  C-D,  D-A)  of  the  car¬ 
bon  plate  are  all  displayed.  Results  show  that  there  exists  an 
appreciable  discrepancy  between  the  simulated  experimental 
data  and  the  numerical  data  of  best  fit.  This  implies  that  the 
convergence  of  temperature  prediction  process  is  relatively 
poorer  when  one  uses  the  temperature  data  on  the  edge  surface 
of  the  carbon  plate.  However,  considering  the  experimental 
cost  reduction,  one  may  still  prefer  using  the  temperature  data 
on  edge  surfaces  of  the  carbon  plate  rather  than  using  those 
measured  on  the  end  plate  surface. 

Effects  of  the  number  of  measurement  points  on  the  edge 
surfaces  of  the  carbon  plate  on  temperature  prediction  are 
also  investigated.  The  predicted  temperature  profiles  along 
y  =  0.105  m  with  different  numbers  of  measurement  points 
for  case  3  at  o  -  0.1  are  shown  in  Fig.  9.  The  number  of 
measurement  points  is  represented  by  (c  x  d)  x  4  where  c 
x  d  represents  the  number  of  measurement  points  on  each 
edge  surface  and  4  accounts  for  the  four  edge  surfaces  with 
a  carbon  plate.  In  this  figure,  the  number  of  measurement 
points  is  assigned  to  be  (21  x  5)  x  4,  (6  x  3)  x  4,  and 
(3  x  1)  x  4,  which  are  also  plotted  in  the  figure.  Again, 
it  is  found  that  the  accuracy  of  temperature  prediction  of  the 
internal  temperature  of  the  fuel  cell  is  significantly  influenced 
by  the  number  of  temperature  measurement  points  on  the 
edge  surfaces  of  the  carbon  plate. 

It  is  important  to  mention  that  the  present  inverse  approach 
is  not  limited  to  the  temperature  measurement  application. 
Distribution  of  the  current  density  and  gas  concentrations 
inside  the  fuel  cell  may  also  be  determined  provided  that  a 
powerful  direct  problem  solver  is  available. 

4.  Concluding  remarks 

The  present  study  is  concerned  with  temperature  predic¬ 
tions  of  the  internal  temperature  of  a  fuel  cell.  The  temper- 
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ature  distribution  at  the  ME  A/carbon  plate  interface  is  pre¬ 
dicted  based  on  the  measured  temperature  data  on  the  outer 
surfaces  by  an  undestructive  inverse  method.  To  demonstrate 
the  validity  of  this  approach,  three  different  cases  (cases  1-3) 
with  different  forms  of  the  exact  solution  are  specified  to  be 
predicted,  and  the  temperature  measurement  data  on  the  outer 
surfaces  are  generated  by  adding  a  perturbation  to  the  exact 
temperature  distribution  on  the  outer  surfaces.  Effects  of  the 
arrangement  and  positions  of  the  temperature  measurement 
points,  the  experimental  uncertainty  ( o ),  and  the  number  of 
measurement  points  on  the  internal  temperature  prediction 
have  been  investigated.  Meanwhile,  the  performance  of  the 
temperature  prediction  based  on  the  edge  temperature  data 
of  the  carbon  plate  is  compared  with  that  based  on  the  outer 
surface  temperature  data  of  the  end  plate.  In  addition,  an  er¬ 
ror  norm  (s)  and  a  specific  objective  function  (j)  have  been 
defined  in  order  to  quantify  the  accuracy  and  convergence  of 
the  approach. 

Four  possible  arrangements  of  the  measurement  points  are 
installed  on  the  outer  surface  of  the  end  plate,  and  the  pre¬ 
dicted  solutions  of  the  third  and  the  fourth  measurement  point 
arrangements  appear  to  be  identical  to  the  exact  solution. 
However,  the  third  measurement  point  arrangement  results 
in  a  greater  error  in  prediction  than  the  fourth  arrangement. 

It  is  found  that  at  a  —  0,  an  accurate  temperature  prediction 
can  be  obtained  for  case  1 .  However,  when  a  higher  uncer¬ 
tainty  of  a  =  1.0  is  introduced  into  the  experimental  tem¬ 
perature  data,  the  deviation  between  the  predicted  and  the 
exact  temperature  data  is  greatly  increased.  The  error  norm  e 
reaches  approximately  0. 07 3  °C  at  a  -  1.0. 

For  cases  2  and  3,  even  at  a  -  0,  the  exact  temperature 
curves  cannot  be  fitted  exactly  by  the  polynomical  tempera¬ 
ture  distribution  like  case  1 .  The  results  imply  that  the  magni¬ 
tude  of  error  is  also  dependent  on  the  features  of  temperature 
distribution  at  the  ME  A/carbon  plate  interface. 

In  addition,  the  comparison  between  the  best  fit  numerical 
data  and  the  simulated  experimental  data  for  the  temperature 
distribution  on  the  outer  surface  shows  that  a  best  fit  can  be 
accurately  yielded  for  each  set  of  the  simulated  experimental 
data.  From  the  specific  objective  function  at  various  experi¬ 
mental  uncertainties  for  cases  1,  2,  and  3,  it  can  be  observed 
that  an  increase  in  uncertainty  leads  to  a  higher  j,  which  rep¬ 
resents  a  worse  convergence. 

Meanwhile,  11x11  measurement  points  are  sufficient 
to  provide  an  accurate  prediction  of  the  ME  A/carbon  plate 
interface  temperature.  The  discrepancy  in  the  predicted  tem¬ 
perature  distributions  between  the  exact  and  the  predicted 
curves  is  increased  as  the  number  of  measurement  points  is 
reduced.  As  the  number  of  measurement  points  is  decreased 
from  21  x  21  to  3  x  3,  the  error  norm  e  is  elevated  from  0.01 
to  0.073  °C. 

As  for  the  temperature  prediction  based  on  the  edge  tem¬ 
perature  data  of  the  carbon  plate,  it  is  found  that  the  con¬ 
vergence  of  temperature  prediction  process  is  relatively  poor 


when  one  uses  the  temperature  data  on  the  edge  surfaces  of 
the  carbon  plate,  instead  of  the  temperature  data  on  the  end 
plate  surface,  and  the  accuracy  of  temperature  prediction  of 
the  internal  temperature  of  the  fuel  cell  is  significantly  influ¬ 
enced  by  the  number  of  temperature  measurement  points  on 
the  edge  surfaces  of  the  carbon  plate. 
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